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1 Introduction 



The Balitsky-JIMWLK equation [1-4] are equations governing the small-x QCD evo- 
lution for dense partonic systems. The Balitsky equations are an infinite hierarchy 
of coupled equations expressing the energy dependence of the scattering of high en- 
ergy quarks and gluons (represented by Wilson lines in the fundamental and adjoint 
representations, respectively) on a target. The JIMWLK equation is a functional 
Fokker-Planck equation [4, 5] for the small-x evolution of the target wavefunction 
equivalent to the Balitsky equations. The Kovchegov [6] equation is a simplified 
version of the Balitsky equations where correlations are suppressed, leading to a 
relatively simple nonlinear equation for the elastic scattering amplitude. 

It has recently [7, 8] been realized that the Balitsky-JIMWLK equations miss 
some essential ingredients in satisfying unitarity constraints in a realistic manner. 
While these equations accurately handle the recombination of gluons when the gluon 
occupation number is large they do not properly create the growth of the occupation 
number starting from a dilute system. For that reason they are accurate, at least 
in a limited energy domain, when starting with a dense wavefunction, such as that 
of a big nucleus, but they are not accurate starting from a dilute system such as an 
elementary dipole. 

Iancu and Triantafyllopoulos [7, 8] suggested a new equation which consists of 
the Balitsky hierarchy along with a stochastic term which, in a dipole language, 
corresponds to dipole creation or dipole splitting. In Ref. [9] Mueller, Shoshi and 
Wong cast this equation into an equation for the JIMWLK weight function with the 
addition to the usual JIMWLK Fokker-Planck term being a fourth order functional 
derivative. Later on in the paper, this extension to the JIMWLK equation will be 
referred to as the MSW term. The effect of this stochastic term on the saturation 
momentum and the scattering matrix at asymptotic rapidities has been worked out 
in [10, 11]. Finally, in Ref. [12] both the splitting and recombination terms were 
written in a compact and simple form in the relevant large N c limit for the high- 
energy scattering problem. (Parts of these results were anticipated in Ref. [13].) 

Kovner and Lublinsky [14, 15] have suggested a general duality between the equa- 
tions for low-density and high-density systems. In the above mentioned large N c 
limit this duality is manifest [12] and the equations for dipole splitting are exactly 
those given in Ref. [8, 9]. The general form of the duality suggested in Refs. [14, 15] 
is likely correct and there is an ongoing effort to make this duality more explicit [16]. 

In this paper, we address the relationship between the splitting term discussed in 
Ref. [8,9] and the more general splitting terms given in Ref. [15]. This is partly an 
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elaboration of discussions previously given in Refs. [17-19] in a closely related con- 
text. In particular we emphasize that the splitting term used in Refs. [8, 9] actually 
defines what is meant by the large- A^ limit in a high-energy scattering problem. We 
recall that for the high-energy scattering of two dipoles the leading terms in a s N C Y 
are in fact just the BFKL parts of the evolution. All multi-BFKL evolutions (multi- 
pomeron terms) are, strictly speaking, higher order in N c . However, it is convenient 
to define a high-energy large- N c limit where one keeps 1/A" C 2 terms which are en- 
hanced by a factor exp[(ap — 1)Y] [17-19] where Y is the rapidity of the scattering 
and «p the usual hard pomeron intercept. This is a natural definition since it is 
exactly when a 2 s exp[(«p — 1)Y] is the order of one that unitarity corrections become 
important. In this large- N c limit there is no difference between the dipole splitting 
used in Ref. [8, 9] and the more general splitting given in Refs. [14, 15]. 

While from the point of view of an evolved wavefunction the dominance of the 
"large-Ay dipole splittings is manifest it is not so clear how this occurs as one 
evolves the general vertex (dipole splitting). We examine both generally and explic- 
itly the scattering to two projectile dipoles on a target dipole. At low energies the 
general splitting term and the "large- N c " splitting term are not the same and the dif- 
ference between the two contributions is of the same size as either contribution [15]. 
However, as one evolves to large rapidity the dominant contribution, involving the 
"large- A^c" splitting term, behaves as afexp[2(a:p — 1)Y] while the difference be- 
tween the large- A" c term and the general term behaves as a A s exp[(a:p — 1)Y]. This is 
a remarkable result. It has a counterpart in pomeron splitting language where Braun 
and Vacca [20] observed that in a particular, and natural, scheme for defining the 
triple pomeron vertex [21,22] that vertex dominates over non-triple pomeron split- 
tings. This dominance seems to be exactly the same as what has been described in 
comparing our large- A" c splitting and the correction terms which evolve less rapidly 
in rapidity. 

We have also found a surprising result. If one writes the evolution of the 
JIMWLK weight function as in Eq. (5.2), then we find X to be scheme depen- 
dent while X Wy is uniquely defined. We have arrived at this result in an explicit 
way. We take X to be given by Eq. (5.10), and we find that the terms involving four 
derivatives with respect to p do not agree with the result, Eqs. (2.38) and (2.39), 
of Ref. [15]. However, we do not believe this to be because one of the calculations 
is incorrect, but rather because the definitions of X are different in the two cases, 
Eq. (2.33) in [15] and Eq. (5.10) here. We believe both definitions of X are cor- 
rect and that when a physical quantity, such as the right hand side of Eq. (5.2), 
is calculated agreement will be reached. We have explicitly evaluated all the terms 
involving fourth order derivatives of p on the right hand side of Eq. (5.2), and we 
find the generalized dipole picture of Kovner and Lublinsky[15] emerges. This can 
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be viewed as an explicit check of the more general derivation of evolution of the 
generalized dipole picture of the weight function W Y given in Sec. 3. (This result 
has also recently been found in Ref. [23].) 

The outline of the paper is as follows: In Sec. 2 we briefly review JIMWLK 
evolution and the dual evolution discussed in Refs. [14,15]. In Sec. 3 we show 
that the dual evolution follows that of a generalized dipole picture as previously 
suggested by Kovner and Lublinsky [15]. In Sec. 4 we show that the splitting 
terms not included in the large- N c picture evolve much more slowly in energy than 
those included in the high-energy large-iV c approximation. In Sec. 5 we explicitly 
evaluate terms through four derivatives in p and discuss the scheme dependence of 
the evolution kernel. 



2 The evolution equation of the color glass con- 
densate 

2.1 The Balitsky- JIMWLK equation 

For a left moving projectile with lightcone time x~ which experiences at the time of 
the collision the effect of the color fields a a (x~,x) originated from the target, the 
Balitsky- JIMWLK equation is 

w w ^ - \ L fo^y y)[a] ^k WAa] • (21) 

where Y = lnl/x is the rapidity of the small- £ gluons, J = J d 2 x (x, y denote 
transverse coordinates), a a is the color gauge field radiated by the color sources 
in the target, and the functional WV[ck] is the weight function for a given field 
configuration. The functional derivatives S/5a a (x) here are meant to be derivatives 
at the largest x~ or equivalently they stand for 

6 = lim >• . • ( 2 - 2 ) 



Sa a (x) x--^oo5a a (x ,x 
The kernel r] ab (x, y)[a] is a non-linear functional of a, 

r) ab (x,y)[a] = - [ IC(x, y, z) (1 + V±V V - V^V Z - V}V y ) ab , (2.3) 

as it depends on the Wilson lines in the adjoint representation V and V\ 

V£[a] = Pexp^ J dx-T a a a (x-,x)^ (2.4) 
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with (T a ) bc = —if abc , where P denotes the path-ordering along the light-cone coor- 
dinate x~ and 

y^) = ihi r " z w (y " % ■ ( 2 - 5 ) 

(2tt) 2 (x - zf{z - yf 



2.2 The dual equation to Balitsky-JIMWLK 

Given the JIMWLK equation, the dual equation can be obtained from JIMWLK 
via the duality transformations [12, 14-16] 

S ^ 
x~ -<=>- x + , - — — — r •<=>- ip a (x + , x) , a a (x~, x) •<=>- —%- 



Sa a (x~,x) Sp a (x + ,x) 

(2.6) 

This transformation is a consequence of the Lorentz boost invariance and the arbi- 
trariness in the labelling of which of the incoming hadrons is the target and which 
is the projectile. The resulting dual equation to JIMWLK is 



_d_ 

dY 



W Y [p] = -\ [ p a (x)r ] ab (x,y)[-t5/5p]p b (y)Wy[p] (2.7) 

A Jx,y 



where p a (x) is now the two-dimensional charge density of the target at the largest 
x + , so similarly to Eq. (2.2) the following short notation is used 

p a (x) = lim p a (x + ,x)= lim / dx~ p a (x + , x~, x) . (2.8) 

x+— >oo x+-^oo J 

The new r] ab [— i5/5p] is given by 

r) ab (x,y)[-i6/6p} = - [ ^x,y,z)(l + UiU y -UlU z -UlU y r b , (2.9) 

Jz 



with 



Ul = V*[-i6/6p] = Pexp^ J dx + T a 5/5p a (x + ,x^J . (2.10) 

Here P indicates, as usual, an ordering of the integral where terms having larger 
values of x + come further to the left. Whereas in the JIMWLK equation the em- 
phasis is on the functional variable a(x~, x), the field experienced by the projectile, 
the dual equation has the emphasis on the two-dimensional density, p(x + , x), of the 
target. 
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3 The color dipole version of the dual equation to 
JIMWLK 



In order to arrive at a color dipole version of Eq. (2.7) upon which all our further 
discussions will be based, we examine two cases of the Wy for a single dipole and 
for two dipoles respectively to all orders in derivatives in p. At the end we make a 
brief remark on the case of an arbitrary number of dipoles. 

3.1 The case of a single dipole 

We denote a single dipole as 

W<?\x,y) = ^-HUlUy) 5[p] = R(x,y)5[p] . (3.1) 

U\. is given by an identical expression as \J\, in Eq. (2.10) but in the fundamental 
representation of SU(N C ). The superscript (0) in Eq. (3.1) stands for the fact that 
one begins with an unevolved dipole, with F-evolution being given by Eq. (3.2). 
Thus at first order in the evolution 

-^W$\x,y) = ~\j u p a (u)r ] a \u,v)[- l 5/5p]p b (v)W^\x,y) . (3.2) 

The action of p a (v) on Wy\x, y) is 

p\v) Wf(x, y) = |- tr{UlU y t b ) (5 2 (v - y) - 5 2 (v - x)) 5[p] (3.3) 

and 

p a {u)p\v) wP(x, V) = Y C {- ^(UlJJ v th a )5 2 {u -x) + tr(UiU y t a t b )5 2 (u - „)) 

x (5 2 (v -y)- 5 2 {v - x)) 5[p] . (3.4) 

To use Eq. (3.3) and Eq. (3.4) one has first to bring p a (u) to the right of rj ab (u, v). 
The action of p a {u) on (U v ) cb is 

p a (u)(U v y b = (U v y b p a (u)+g(U v T a y b 5 2 (u-v) . (3.5) 

Using Eq. (3.3) and Eq. (3.4) in Eq. (3.2) one obtains 

-^w^ix, y)= ^fj z v,*){i + ulUy - uiu z - uliiyY b 
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xtr(U x U y t b t a )5[p] 
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c 

2 



2nN c j z 



K(y, y, z) {2 - tj y \J z - UtU y } ab tr{U x U y t b t a ) S\p] 
[ ()C(x,x,z){UlU x T a } ab -lC(y,y,z){UiU y T a } a < 



xtr(UtU y t b )6\p] . 

(3.6) 

To proceed further we note that U x can be expressed entirely in terms of U x 

{U x } ab = 2ti{Uit a U x t b ) (3.7) 
and one can use the Fierz identity 

= ^ (M* " ;^M«) (3-8) 

so the product of C/'s can be rewritten as 

{UlU y } ab = {Ul} ac {Uy} cb = 2Hr(Ult c U x t a )tr(U y t c U y t b ) = 2tr(U y U x t a UlU y t b ) . 

(3.9) 

Now repeated use of Eq. (3.8) allows us to work out each term in Eq. (3.6). For 
example 

{U x U y Y b tr{U x U y tH a ) = \{N C - ^) tr{U x U y ) , (3.10) 

{u x u z y b t?{u x u y t b t a ) = {UtU v Y b tr{U x U v t b t a ) 

= l -ti{U x U z )tT{UlU y ) -±-tr(UlU y ) , (3.11) 

{UtU x } ab tT(UtU y t b t a ) = \tr(UiU x )tT(UlU z UiU y ) - JLtr(E/ltf w ) . (3.12) 
In the last line of Eq. (3.6) there are terms with 

{UlU x T a } ab = 2ti{U x U z t a UlU x t c ) {T a } cb = -2tr{U x U z t a UlU x {tH a - t a t b )) 

= - tr(UtU z )tr(UtU x t b ) + ti{UtU x )ti{U x U z t b ) . (3.13) 

This contracting with the remaining trace factor gives 

{UlU x T a y\r{U x U y t b ) = ~{tr(UtU z MUlU y ) - trfE/JOK^.^)) . 

(3.14) 



Substituting everything into Eq. (3.6) all terms accompanied by an explicit l/N c 
factor such as the second term in Eq. (3.10), Eq. (3.11) and Eq. (3.12) cancel within 
each {} ab and terms with a trace of four C/'s also cancel leaving 

^ W$\x, y)=-J^- J (-2K(x, y, z) + JC(x, x, z) + K(y, y, *)) 

9 2 N t 
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M{x, y, z) {R(x, y) - R(z, y)R(x, z)} 5[p] 

(3.15) 



where 



is the kernel of the dipole model. Eq. (3.15) is valid for any values of N c and is the 
equation reminescent of the dipole form of the Balitsky-Kovchegov equation that 
most of the following discussions will be based. 



3.2 The case of two dipoles and beyond 

Now starting with two dipoles, that part of Wy is 

W<?\x,w,y) = ^(UlXJj^(UlU y )5[p] = R(x,w) R(w,y) S[p] . (3.17) 

Once again we evolve this by using Eq. (2.7) or Eq. (3.2) and eliminating first the p's 
by acting them on Wy\ As before there are two possibilities. The first possibility 
is only one p acts on Wy^ 

p\v)wP{x,w,y) = R(x,w) (p b {v)W^\w,y))+R{w,y) (p\v)W$\x, w)) 

(3.18) 

while the other on rj ab , where the (pWy^)^ on the right hand side are given by 
Eq. (3.3), and the second is both p's act on W Y °^ 

p a (u)p b (v)W Y d) (x,w,y)= R(x,w) (p a (u)p b {v)W$\w,y)) 

+ R(w,y) (p a (u)p h {v)Wy"\x,w)) 

+ (p\u)R(x, w)) (p\v)R{w, y)) 5[p] 

+ (p b (v)R(x, «;)) (p a (u)R(w, y)) 5[p] . (3.19) 
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The two {ppW^ys are given by Eq. (3.4) and the (p-R)'s are essentially the same 
as (pWy^) in Eq. (3.3) but without the explicit 5[p] factor. These are of course 
related by 

(pR(x,y))5\p] = (pW$\x,y)) . (3.20) 

Based on what we learned in Sec. 3.1 and a careful examination of the form of 
Eq. (3.18) and Eq. (3.19), it is easy to see that 

r\ r\ r\ 

— W?(x, w, y) = R(w, y) — W?\x, w) + R(x, w) — w£\w, y) 

- l -f v a \u,v)[-i5/5p]^ a (u)R(x,w))^ b (v)R(w,y))5[p] 

- l -J v ab (u,v)[-i5/5p](p b (v)R(x,w)) (p a {u)R{w,y)) 5[p] . 

(3.21) 

The expressions in the last two lines can be worked out. The first of these is 
- l - J v ab (u, v) [-t5/5p] (p a {u)R{x, w)) (p\v)R{w, y)) 5[p] 

= ~ j, K( > W > y > Z ^ 1 + UlUy - UlU Z - UlUyY b 

+ K(x, w, z){\ + UlU w - UlU z - UlU w } ab 

- lC{w, w, z){2 - UlU z - UlU w } ab 

- JC(x, y,z){l + UtUy - U ] X U Z - UlU y } ab ) 

x tr(UlU w t a ) tr(UlU y t b ) 5[p] . (3.22) 

Again with the use of the Fierz identity Eq. (3.8), one can work out the terms. For 
example 

{UlUyYMuiU w ntv{UlU y t b ) = \tv{UlU y ) - ^tr(UtU w )tr(UlU y ) , (3.23) 

{UlU z } ab tT(UlU w t a )tT(Ui,U y t b ) = hi(uiu z ulu y ulu w ) 

' tr(UlU w )tr(UlU y ) , (3.24) 



2N C 



{UlU w } a hi(UlU w t a )tr(UlU y t b ) = l -tr{UlXJ w UlU y UlU z ) 



' tr{UlU w )tr{UlU y ) . (3.25) 



2iV c 



g 2 



In fact if one ignores the /C factors, all terms in Eq. (3.22) with {UW z } ab give the 
same expression on the right hand side as Eq. (3.24), all terms with {U z U} ab give 
Eq. (3.25) and lastly all terms without either U z or U\ give Eq. (3.23). Substituting 
all these back into Eq. (3.22), terms with a 1/N C factor cancel among themselves 
within each {} ab and one is left with 

- l -J r] ab (u,v)[-i5/5p}[p a (u)R(x,w))[p b (v)R(w,y))5\p} 

J (fc(w, y, z) + K,(x, w. z) - K,(w, w. z) - K.(x, y, z)j 

x (2tr(UlU y ) - tr(UlU z UlU y UtU w ) - tr(E£tf w E^^)) &\p\ 

(3.26) 

Incidentally swapping p a and p b in the above equation gives exactly the same result 
so the last two lines of Eq. (3.21) are identical. Finally Eq. (3.21) becomes 

d 

— W^\x,w,y) 
= R(w, y) — W£\x, w) + R(x, w) — W$\w, y) 
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J (fc(w, y, z) + fC{x, w. z) - K{w, w, z) - fC{x, y, z) j 
x (2tr(f4ig - tr(UlU z UlU y UtU w ) - ti{UlU w UlU y UlU z ) ) 5[p] 



R(w, y) — W$\x, w) + R(x, w) — W?\w, y) 
9 2 



2ttN c j _ 



(M(x, y, z) - M{w, y, z) - M(x, w, zfj 

x (2R(x,y) - ^[tr(UiU z UlU y UlU w )+tr(UlU w UlU y UlU z )^ 5[p] 

(3.27) 

Apart from the gauge field A^s in the C/'s and W's have been replaced by S/5p here, 
the presence of S[p] and a prefactor, this last expression is the same equation as Eq. 
(121) of the first paper in Ref. [1]. Using Eq. (3.15) one can see that the terms in 
the second line is down by a factor of l/iV c 2 in comparison to those in the first line. 
Therefore in the large iV c limit, the evolution equation for two dipoles simplifies to 
a "product rule" like form 

r\ r\ r\ 

— (x,w,y)~R(w,y) — (x,w) + R(x, w) — (w,y) . (3.28) 

10 



Because there are only two p's in Eq. (2.7) and they can act on a maximum of two 
different i?'s, this product rule holds for the evolution of any arbitrary number of 
dipoles in the large N c limit. One can conclude that in this limit each dipole evolves 
independently of each other. 

4 Dipole model versus dual evolution in the low 
gluon density region 

In this section we consider the scattering of two elementary dipoles off an evolved 
target. The latter is considered to be dilute (low gluon density) so that gluon 
mergings can be neglected. We construct the wavefunction of an evolved target 
starting with an elementary dipole by using the dipole model [24, 25] and the dual 
evolution equation [14-16]. With the evolution operator recently worked out by 
Mueller, Shoshi and Wong (MSW) [9] to extend the JIMWLK equation to the low 
gluon density region, we obtain a scattering amplitude which in the dipole model is 
dominated by the scattering of the two projectile dipoles off two different dipoles in 
the evolved target. Kovner and Lublinsky [15] have shown that the dual evolution 
equation gives additional terms, missed by the MSW evolution operator, which 
allow for the two projectile dipoles also to scatter off a single dipole in the evolved 
target. In this section we show that the multiple interaction of the target dipole is 
suppressed as compared to the single interaction of different target dipoles at high 
rapidities. The two contributions are comparable only at very low rapidities where 
the target contains only few gluons. At high rapidities with a large number of gluons 
in the target the effective operator of Ref. [9] gives correctly the leading contribution 
to the scattering amplitude. 

4.1 Low gluon density evolution in the dipole model 

Let us consider the evolution of an unevolved single dipole in the weak field limit 



Here the subscript (2) stands for the second order derivative in p and the superscript 
(0) denotes the unevolved state of the dipole that one begins with. The 8 b x here is 
given by 



(4.1) 




(4.2) 
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The first step in the evolution in the extended JIMWLK equation as given by the 
dipole model reads [9] 

W? = (X 2 + XT W ) < ( } 2) , (4.3) 

where the operator X2 gives the BFKL evolution and X4 ISW splits one scattering 
dipole into two scattering dipoles at every evolution step. In the JIMWLK language, 
in terms of fields, these have been calculated in Ref. [9]. Using the relation between 
the fields and charge densities — a a (x + , x) = p a (x + , x), they can be expressed in 
terms of charge densities 

*' = 2^^ L VZ M{U ' V ' Z) K^-^-'ft] (T c T d ) ab p a {u)p\v) 

(4-4) 



and 

X 



1 a 2 N 

msw _ x y jv c 



/ M(u,v,z) [{5 c u -8%f] [(5 d z -5 d v y]p a (u)p%v) 



_d_ 

dY 



16g 2 N* 2tt 

(4.5) 

Acting with these operators on the Wy^(x,y) in Eq. (4.3), one obtains 
W$\x,y) = ^ jM(x,y,z) (4.6) 



x 



z) + R^(z, y) - R%{x, y) + R$(x, z)R { § ) { 



where the first three terms on the right-hand side (r.h.s.) represent the BFKL 
evolution of the original dipole (x, y) while the fourth term describes the splitting 
of the original dipole into two new dipoles (x, z) and (z, y) which simultaneously 
scatter off the projectile. This first step in the evolution of an unevolved single 
dipole is shown in Fig. 1. 

Let us now consider the second step of the evolution. A graphical representation 
of the result is shown in Fig. 2. This is simply obtained by evolving one step further 
the graphs on the r.h.s. of Fig. 1 according to the dipole model in the extended 
evolution, i.e., by acting on them with X2 and X4 SW as in the first step. In Fig. 2 
the first five lines come from the BFKL evolution of the four graphs in Fig. 1, 

(< } 2) + w^n with 

Wf ( y w (*, V, z) = B$(x, z)R$(z, y) 5[p] . (4.7) 

The last line results by acting with X4 ISW on the first three graphs in Fig. 1, 
x msw W W Tlie su b S cript [ ] 4 in Fig. 2 and Eq. (4.8) indicates that we have 
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Figure 1: One step of evolution of an unevolved single dipole in the dipole model: 
The first three graphs on the r.h.s. represent the BFKL evolution of the original 
dipole and the last graph the splitting of the original dipole into two new dipoles 
which simultaneously scatter off the projectile dipoles. 



dropped in the wavefunction terms having more than four orders in p-derivatives, 
like Xf sw Wyj^ SVf , since for our purpose it is enough to consider the scattering of 
the evolved target on only two projectile dipoles. 

In formula the second step of the evolution is 

w K 2) ] 4 = X2 w ^ + X2 w ^ sw + xfSW w ^ • (4 - 8) 

It is easy to check that this evolution equation holds for any step of the evolution. 
Summing over all dipoles numbers (superscript), the evolved wavefunction which 
contains up to four p-derivatives reads 

^{WyI = ^W y{2) + ^W^ (4.9) 

with 

d 

—W Y ( 2 ) = X 2 W Yi 2) , (4.10) 

In Eq. (4.11) Xf sw W Y (2) converts one into two simultaneously interacting target 
dipoles at any step of evolution. Both dipoles evolve in the subsequent steps of the 
evolution according to the BFKL equation, given by X2 Wy ST, before they scatter 
on the two projectile dipoles. If the two dipoles are created in the first step of the 
evolution (fourth graph on the rhs of Fig. 1), then the subsequent BFKL evolution 
of both dipoles gives for the rapidity dependence of the wavefunction 

with ap — l = (4a s N c In 2) / ir. When considering the scattering of the evolved target 
dipole off the two projectile dipoles, (4.12) gives for the T-matrix 

T 2Pom (F) ~ at e 2{ap ~ 1)Y . (4.13) 
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Figure 2: One step of evolution of the graphs on the r.h.s. of Fig. 1 in the dipole 
model. More than four gluon interactions with the projectile dipoles are excluded. 



In the Pomeron language the interaction of the evolved target dipole with the two 
projectile dipoles occurs via the double Pomeron exchange as shown in Fig. 3(a). 
However, if the two target dipoles are generated after several previous steps of 
BFKL evolution of a single dipole (last line in Fig. 2), then the rapidity depen- 
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dence is partly that of a single Pomeron and partly that of a double Pomeron ex- 
change, e^ ap ~ 1 ^ y e 2 ^ ap ~ 1 ^ Y ~ y \ which is suppressed as compared with the pure double 
Pomeron exchange. The later case is shown in Fig. 3(b). 




Figure 3: An evolved target scattering off two elementary dipoles in the Pomeron 
language: (a) double Pomeron interaction, (b) a partly evolved Pomeron splitts into 
two Pomerons which after some evolution scatter off the two dipoles, (c) a fully 
evolved Pomeron scatters at the end of its evolution off the two dipoles. 



4.2 Low gluon density evolution in the dual equation 



The first step of evolution of an unevolved single dipole as given by the dual equation 
Eq. (2.7) up to the fourth order in p-derivatives (see Eq. (3.2) and Eq. (3.15)) is 



_d_ 

dY 



(1) dual 
Y 



d_ 

9 2 N ( 



Y(2) 



w (1) +w {1) 

VV Y(3) + VV Y (4) 



J4 



(4.14) 



2tt 

x 



M(x,y,z) 



x,y,z 



R^(x,z) + R^(z,y)-R^(x,y) 
+ R^(x,z) + R^(z,y)-R ( § ) (x,y) 
+ R^(x,z)R < § ) (z,y) 
+ R$(x,z) + i$j(z,y) - R<g>(x,y)]6\p] 



with denoting the n-th order p-derivative term in the expansion of W?> m 

Eq. (3.1). Fig. 4 shows Eq. (4.14) in graphical form. Note that the dual evolution 



,(o) 
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allows for a single target dipole to interact also via three and four gluon exchange 
while in the dipole model a single target dipole interacts only via two gluon exchange. 

The second step of the evolution follows simply by applying the same rule as in 
the first step for the created dipoles inside the evolved target dipole. The resulting 
graphs are shown in Fig. 5. As compared with the dipole model, the extra graphs 
in the dual evolution are those with three and four gluons attached at the end of 
the BFKL evolution to a single target dipole. Further steps of evolution will only 
add some more dipoles according to BFKL but will not split the three and four 
gluons and attach them to different target dipoles. Thus, the rapidity dependence 
of the evolved target wavefunction where three or four gluons are attached to a 
single target dipole, Wy, is that of a single BFKL Pomeron 

W Y(l) ~ e {ap -^ Y W Y=0{t) (4.15) 

where % = 3 or 4 stands for three or four gluons respectively. This leads to the 
following rapidity dependence of the T-matrix for the scattering of the two projectile 
dipoles via two-gluon exchange with a single dipole in the evolved target 

T 1M (y) ~ a 4 s e {ap ~ 1)Y . (4.16) 

In the Pomeron language this result is shown in Fig. 3(c) where a single Pomeron 
interacts at the end of its evolution via four gluon exchange with two dipoles. Note 
that the three gluon exchange part Wy(z) will drop out when considering the inter- 
actions with two projectile dipoles. 

The wavefunction of an evolved target dipole at the fourth order of p-derivatives 
as given by the dual equation can be separated into two parts in the form 

=W?l$ + W Y{A) , (4.17) 

where the first part is in the dipole model and the second is not. Eq. (4.17) tells 
us that at high rapidities Wy ^ gives the dominant contribution to the scattering 

matrix since Wy S ^ > Wy (4) (see Eq. (4.12) and Eq. (4.15)). This is the same result 
that was recognized already in Ref. [20]: the double Pomeron interaction (Fig. 3(a) 
and Eq. (4.13)) dominates over the single Pomeron which interacts at the end of its 
evolution with the two external dipoles (Fig. 3(c) and Eq. (4.16)). 

In the dipole model [24, 25] our main result reflects the fact of dipole counting in 
the scattering of two external dipoles with the evolved target: when both external 
dipoles can scatter off different target dipoles, each of them participates with the 
factor a 2 n(Y), thus, giving T ~ [a 2 n(Y)} 2 , with the dipole number density n(Y) ~ 
exp [{dtp — l)Y]. On the other hand, when both external dipoles are forced to scatter 
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Figure 4: One step of evolution of an unevolved single dipole according to the dual 
equation up to the fourth order in p-derivatives. 

off a single target dipole, the one scattering first gives the factor a 2 s n(Y) while the 
second one has only one possibility a 2 x 1, leading to T ~ \a* n(Y)]. Note that only 
at very low rapidity where n(Y) is of order one, the two cases give similar results. 



5 Scheme dependence of the evolution kernel; ex- 
plicit calculations 

In Ref. [15] the evolution in rapidity of the target wave function was written in 
the Hamiltonian form. Two versions of this Hamiltonian were given explicitly in 
terms of derivatives of p up to the fourth derivative: one for the general case and 
the other for when one is dealing exclusively with dipoles. The problem is that 
the derivative form of these two Hamiltonians did not appear to be equivalent to 
each other. On top of that the dipole version of this Hamiltonian was emphatically 
pointed out to differ from the MSW term. The reason for the latter was partly 
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d 

dY 



graphs on the 
rhs of Fig. 4 



graphs from the dipole model 
on the rhs of Fig. 2 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ 



+ = =, 



Figure 5: One step of evolution of the graphs on the r.h.s. of Fig. 4 according to 
the dual evolution up to the fourth order in p-derivatives. 
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explained in Sec. 4. There remains the disagreement of the general derivative form 
and the dipole derivative form of this Hamiltonian in [15]. It is the purpose of this 
section to show that one should not compare these Hamiltonians among themselves 
as well as comparing them with the MSW term directly because the specific form of 
the Hamiltonian itself is not unique hence the scheme dependence (see below). For 
this reason it only makes sense to compare physical quantities such as the action 
of these Hamiltonians on the target wave function. Those are the quantities that 
should be in agreement independent of the form of the Hamiltonian. 

From the discussion in Sec. 3 it is clear that our general form of Wy, evolved 
from a parent dipole with transverse coordinate (x,y), is 

oo „ N „ 

W Y (x, y) = J2 dT NPN(Y) n / Xi) 5[p] (5.1) 

N=l i = \JXi-l,Xi 

where R(x,y) is as given in Eq. (3.1) while Pn(Y) gives the probability to find a 
system of N "dipoles" at transverse coordinates (xq,Xi), (xi,x 2 ), ... (xn-i,Xn) 
with the first coordinate of the pair (cCj_i, xi) denotes a "quark" and the second an 
"antiquark" with x = x and x N = y, matching onto the parent dipole. One can 
write the evolution of W Y as 

^=XWy. (5.2) 

In Eq. (3.2) we have taken 

X = -\ [ p a (x)r ] ab (x,y)[- i 6/6p]p\y) (5.3) 
J x,y 

with i] given in Eq. (2.9). While this is a natural choice for X it is not unique. Indeed 
the choice made in Eq. (2.33) of [15] is not the same as we have taken although it 
is an equally good choice. The product XWy in Eq. (5.2) is unique, however there 
are different X's which acting on Wy, give the same result. 

The ambiguity in X is a sort of scheme dependence and can be illustrated by 
considering the action of p a (u + ,u) on U\.. It is then straightforward to verify that 

p a (u + ,u)Ul6[p] = {U u {u\u + )Y h p\u\u)Ul5[p\ (5.4) 

where 

U u (u' + ,u + ) =Pexp(-g J dx + T a 5/5p a (x + ,u)^ (5.5) 

and where P indicates that larger x + -components of the integral are put further to 
the right. However it is not true that 

p a (u + , u) = U« b (u' + , u+) p b (u' + , u) . (5.6) 
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The fact that Eq. (5.6) is not generally correct can easily be checked by considering 
the commutator of the operator, 8 X , defined earlier in Eq. (4.2) (we repeat the 
definition here for convenience) 



on the left and right hand sides of Eq. (5.6). Clearly 

[6 x ,p a (u + ,u)]=g5 ac 6 2 (x-u) (5.8) 

while 

[£, U° b (u' + : u + )p\u + , u)} = g 5\x - u) £>>'+, u + ) (5.9) 
and the right hand sides of Eq. (5.8) and Eq. (5.9) are not the same. 

To see this scheme dependence more explicitly we now evaluate X through terms 
including four factors of S/Sp, the same level of accuracy as kept in [9] and in [15]. 
From Eq. (2.9), Eq. (5.2) and Eq. (5.3) 

X = ~f K(x,y,z)p a (x)M ab p b (y) (5.10) 

A7T Jx,y,z 

where 

M ab = {1 + U x U y - U X U Z - U z U y } ab . (5.11) 

To simplify the discussion a little we suppose that X Wy will be used to evaluate the 
scattering of two separate dipoles on the target. In that case M ab can be expressed 
in terms of the operators, 8 X , defined in Eq. (4.2), while in the general case the 
x + -dependence of derivatives with respect to p must be kept in M ab . Then it is 
straightforward to get 

M (2) = _ {TCTd}ab _ 5%) {s , _ ^ (5 12) 

= - \{T c T d T*Y b \5 c x (5 d x - 5 d y )5 y - 5 x (5 d x - 5 d z )5 z - 5 c z (5 d z - 5 d y )5 y ] (5.13) 



^ ^rj^Crj-idrj-ierj-i f ^ab 



$ c J d J e Jl ~ 25 c J d z 5 e J f z - 2S c J d J e z 5 f y + 35 c Jt5 e Jl 



J x^ z^ z" z 



'X X z z 



+ mX^y - WJiStft - 25 z 5 d 5 e y 5f y + 25 x 8 d J x 8 y 



+ 25 c x 5 d y 5 y 5l - 35 x 5 d J y 5 y 



(5.14) 



A casual glance at Eq. (5.14) reveals that this is not the same expression as that 
given by Eq. (2.38) in [15]. In particular Eq. (5.14) has terms which have no 5 Z 
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factors while the formula in [15] has at least one factor of 5 Z in each of its terms. 
As we have no reason to expect agreement when comparing comparable terms in X 
itself, the lack of an agreement is not surprising. 

We now turn to the product XW Y which is not scheme dependent, and we shall 
see how the factors in Eq. (5.12), Eq. (5.13) and Eq. (5.14) organize thermselves 
into a simple expression when the appropriate factors in W Y are included. To make 
the calculation easier we take Wy as 

= R(x,y)5[p]. (5.15) 

Because of the structure of Eq. (5.1), this is no essential limitation. Then letting 
JC{x,y,z) -> -\M(x,y,z), 

X Wf } = - J M (x, y, z) M ab N ab 5[p] (5.16) 

and 

N ab = tr(t a UlU y t b ) . (5.17) 

It is straightforward to evaluate N ab at the lowest few orders in the number of 
p- derivatives. The results are 

NS = 3S* (5-18) 

N^^tr(t b t a t c )(5 c x -5 y ) (5.19) 
1 
2 1 



= ^(t b t a t c t d ) M + *yK - mi) (5.20) 



where, as for M ab , the superscript denotes the order of the derivatives in p. In the 
product of M ab N ab appearing in Eq. (5.16) the term M (2) iV (0) gives the usual BFKL 
evolution 

M m N a b (0) = _ _ (5 



M W N ab (o) vanis hes and M^N^ gives the BFKL evolution of the odderon 

M$N ab W = d abc (s a J b y (5 c x - 5 y ) - S x S b z (S x - 5%) - 5 a z 5 b y (5 z - 5 y )) (5.22) 
while the combination 

M^N^ + M%>N% + M%>N% (5.23) 
gives the four derivative terms which are our main concern here. 
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Using, at large N c 



tr(t b t a t c t d ) 


= ^ (8ab$cd + S bd 5 ac ) 


(5.24) 


(T e T f ) ab tr(t b t a t c t d ) 


1 1 

— -^OcdOef + gd ce % 


(5.25) 


tr(T c T d T e T f ) 


= S c dS e f + 5 c f5de + ^ce^df 


(5.26) 


(T c T d T e ) ab tr(t b t a t f ) 


= - ^tr {T c T d T e T s ) 


(5.27) 



we find 



M W N a b (o) = _ -{5 z .5 z f + -5 Z .5 Z (S Z -S X + S X .S V ) - -5 Z .5 Z (S X -S X + My) 
3 5 

- Y^((5 X .5 Z f + (5y5 Z ) 2 ) + — (S X .S X S X .5 Z + 5y.5y5y.5 Z ) 

5 13 

- 77T^3;-5y(^-5cc + V^y) + q 4c-4c#y -<5y + 777 (^c^y) 2 , (5.28) 

z4 o lb 

(1) = ^-^(M* + My " 2 My) + l(S x .5J y .5 z + 5y.5y4.5 z ) 

3 5 
+ —(5 X .5 Z - 5 y .5 z ) 2 - — (5 X .5 X 5 X .5 Z + 5 y .5 y 5y.5 z ) 

3 5 
+ —5 x .5 y (5 x .5 z + 5 y .5 z ) + —5 x .5y(5 x . 5 X + 5 y .5y) 

1 3 

- ^cc-^y.tfy - -(5 X . 8 y ) 2 , (5.29) 



and 



M (2) Nab (2) = _ + 5yMy _ 2Sx .5 y ) _ -L(5 X .5 Z - 5 y .5 z f 

3 3 
+ —(S X -5 X 5 X .5 Z + 5y.5y5y.5 z ) - — 5 x .5y(5 x .5 z + 5 y .5 z ) 

3 11 
- j^S x .5 y (5 x .5 x + 5y.5 y ) + -5 x .5 x 5y.5 y + -(5 x .5 y ) 2 . (5.30) 

Putting everything together, the result is 

M (4) N ab (0) , M (3) N ab (1) , M (2) ^ab (2) 
ao 1 ao 1 at) 

= ^{10^(5,.^ + 5 Z .5y) " 5(5 Z .5 Z ) 2 - 6^.^(^.^ + 5y.5y) 

- 3((5 X .5 Z ) 2 + (5 y .5 z ) 2 - (5 x .5 y ) 2 ) + A(5 X .5 X 5 X .5 Z + 5 y .5 y 5y.5 z ) 

+ 6(S X .S X Sy.5 Z + Sy.5y5 X .5 Z ) — 125 X .5 Z 5y.5 Z — 4:5 X .5y(5 X .5 X + 5y.5y) 

(5.31) 
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Referring back to Eq. (5.16) one now finds the equality 

9 2 N C 



2tt 



/ M(x, y, z) \r(x, y) - R(x, z)R(z, y)] 5[p\ (5.32) 

Jz 1 J ( 4 ) 



where the subscript (4) in Eq. (5.32) indicates that the term in brackets is to be 
evaluated to fourth order in p-derivatives. We note for completeness that 



(5.33) 



R U^ y^ = W M^OE^ - M a J b J°y + - 8$F V ) (5.34) 



and 



(5.35) 

It is now clear that the fourth order terms in Eq. (5.16) agrees with the general 
expression Eq. (3.15) derived earlier. It is also clear that the simple and elegant 
form embodied in Eq. (3.15) and Eq. (5.32) is not embodied in X alone but is rather 
a property of the product X Wy. 
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